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ABSTRACT 

The details of trajectories of charged particles become increasingly important for proper under- 
standing of processes of formation of radiation in strong and curved magnetic fields. Because of 
damping of the perpendicular component of motion, the particle's pitch angle could be decreased by 
many orders of magnitude leading to the change of the radiation regime - from synchrotron to the 
curvature mode. To explore the character of this transition, we solve numerically the equations of 
motion of a test particle in a dipole magnetic field, and calculate the energy spectrum of magnetic 
bremsstrahlung self-consistently, i.e. without a priori assumptions on the radiation regime. In this 
way we can trace the transitions between the synchrotron and curvature regimes, as well as study the 
third (intermediate or the so-called synchro-curvature) regime. We briefly discuss three interesting 
astrophysical scenarios, the radiation of electrons in the pulsar magnetosphere in the polar cap and 
outer gap models, as well as the radiation of ultrahigh energy protons in the magnetosphere of a 
massive black hole, and demonstrate that in these models the synchrotron, synchro-curvature and 
curvature regimes can be realized with quite different relative contributions to the total emission. 
Subject headings: gamma rays: theory — magnetic fields — radiation mechanisms: nonthermal 



1. INTRODUCTION 

The efficiency of synchrotron radiation depends on the 
pitch angle between the magnetic field and the particle 
velocity. The damping of the perpendicular motion in 
course of radiation reduces the pitch angle. Typically, in 
moderate magnetic fields the pitch angle changes slowly. 
Therefore, for calculations of synchrotron radiation, it is 
sufficient to specify the initial pitch angle distribution of 
particles. 

The situation is different in strong magnetic fields, 
namely when the energy losses become so intensive that 
on fairly short timescales the pitch angle can be reduced 
by orders of magnitude. In this regard, the adequate the- 
oretical treatment of particle trajectories becomes a key 
issue for correct calculations of radiation properties. In a 
curved magnetic field the strong damping of the perpen- 
dicular motion causes transition from synchrotron to the 
curvature radiation regime of radiation. The solutions 
of equations that describe the particle motion with an 
inclusion of energy losses, allow us to trace this transi- 
tion, and thus to calculate the magnetic bremsstrahlung 
without additional assumptions regarding the radiation 
regime. 

In this paper, we study the case of motion of a charged 
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particle in the dipole magnetic field, and calculate self- 
consistently the radiation spectrum taking into account 
the time-evolution of particle's energy, coordinates, and 
direction. We demonstrate that even a small deviation 
of particle's initial direction from the magnetic field line 
may have a strong impact on the character of radiation. 
Despite the fast transition to the final (curvature) regime, 
the particle radiates away the major fraction of its en- 
ergy in the initial (synchrotron) or transitional (synchro- 
curvature) regimes. Consequently, the energy spectrum 
of radiation may differ considerably from the spectrum 
of curvature radiation. 



2. GENERAL COMMENTS 

In a strong magnetic field, the motion of a charged par- 
ticle perpendicular to the field is damped. The amount 
of energy lost during this process depends on the initial 
energy and the pitch angle. In the case of relativistic 
motion along and perpendicular to the magnetic field, 
P|l 3> mc and p±^ ^ mc, particles can lose a large frac- 
tion of their energy even if the initial pitch angle is small. 
Indeed, after the complete damping of the perpendicular 
component of motion in a homo geneous magnetic field, 
the p arallel momentum becomes ( Kelner and Aharonian 
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Then for small initial pitch angles a <C 1 and for the 
relativistic motion perpendicular to the magnetic field, 
we obtain that the final Lorenz factor depends only on 
the initial pitch angle: 
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It is convenient to rewrite this expression in the following 
form: 

(3) 
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where "f± — "f sin a is referred hereafter as perpendicular 
Lorenz factor. We see that at an ultrarelativistic motion 
of particle, even a tiny deflection a = "f±/"f from the 
magnetic field line can cause large energy losses AE ~ 
E{l-lh^)^E. 

The damping rate depends on the strength of the mag- 
netic field. The super-strong magnetic fields surround- 
ing of compact astrophysical objects, such as pulsars and 
black holes, cause very fast damping of the perpendicular 
component of motion, and in this way force the particle 
to move along magnetic field lines. However, since the 
magnetic field lines are curved, the change of the field 
direction results in an appearance of the perpendicular 
velocity. The curvature of magnetic field lines and the 
presence of the perpendicular velocity lead to the cur- 
vature drift (an averaged motion perpendicular to the 
magnetic field due to non-compensated differences in the 
trajectory of the periodic motion arising from changes 
of the magnetic field direction). Thus, after losing most 
its perpendicular motion in a curved magnetic field, the 
particle moves along the drift trajectory gyrating around 
it. 

The radiation spectrum is determined by the curva- 
ture of particle's trajectory which is a superposition of 
the drift and the gyration around it. As long as the real 
trajectory is close to the magnetic field lines, for calcula- 
tions of radiation spectra one can use the curvature of the 
magnetic field lines instead of the trajectory curvature. 
The difference between the magnetic field curvature and 
the drift trajectory curvature can be neglected until the 
curvature drift velocity is small. In case of the magnetic 
field of a long straight wire, the curvature of the drift 
trajectory is 
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where ro is the curvature radius of the magnetic field, 
/3£) is the drif t velocity in the units of the s peed of light 
expressed as (Alfven and Falthammar 1963| ) 
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where = eB/mcy is the gyration frequency, and B is 
the strength of the magnetic field. The gyration itself 
introduces a much larger difference if the velocity per- 
pendicular to the drift traicctory \ is of the order of 
the drift velocity (Kclncr and Aharonian 2012) 
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Note that the particle in the magnetic field of a long 
straight wire can move strictly along the drift trajectory 
without gyration (/3j^ = 0) and therefore with minimum 
possible energy losses. Because of the possibility of such 
motion, an d treating the motion w i th gy ration as a per- 
turbation, Kelner and Aharonian ( 2012 ) has called the 
drift trajectory as a "smooth trajectory" . Having the po- 
tential minimum on the drift trajectory, the particle with 
any initial pitch angle due to energy losses will asymp- 
totically reach the motion strictly along the drift trajec- 
tory. In the case of an arbitrary magnetic field, when 
the curvature is variable, it is not possible to make a def- 
inite statement, except that the particle tries to reach 
the potential minimum according to the local values of 
the magnetic field. If the curvature changes slowly, the 
particle motion could become very close to the drift tra- 
jectory {(3±_ « 0), but, because of gyration, never strictly 
approaches it (/3_l = 0) as in the case of magnetic field 
of an infinitely long straight wire. 

The energy loss rate and the radiation spectrum be- 
have differently from the case of the curvature radiation 
when the velocity component perpendicular to the drift 
trajectory fi±_ is comparable or larger than the drift veloc- 
ity Pd- Because of gyration, the radiation is expected to 
be different from the pure curvature radiation, even when 
initially the particle moves strictly along the magnetic 
field line (/3j_ = /3_d). The nominal curvature radiation 
i s generallv treated as the same synchrotron radiation 
( Landau and Lifshitz| 1975) with a spectral maximum 
at the characteristic frequency corresponding to the 
curvature radius of the magnetic field line instead of the 
curvature radius of the real trajectory. The substitution 
of the real trajectory curvature affects the radiation spec- 
trum, namely shifts the maximum to ~ (1 -I- P^/ I3d)u)^,. 
Thus, the difference from the curvature radiation be- 
comes significant when (5±/Pd ^1- If the pitch angle 
is much larger than the angle between the magnetic field 
line and the drift trajectory, one can neglect it and con- 
sider P± as the velocity perpendicular to magnetic field 
line (although (3±_ is the velocity perpendicular to the 
drift trajectory). It implies the dominance of the ratio 
Pa^/Pd in the expression 1 -|- P±/Pdi i-e. the radiation 
deviates significantly from the curvature radiation when 
the velocity perpendicular to the magnetic field is greater 
than the drift velocity. 

It is interesting to study the possibility of of pitch an- 
gles greater than Pij. This question is specific and its 
answer depends on the acceleration mechanism, but here 
we would like to discuss very general points. The pitch 
angle of the particle accelerated by the electric field in 
the presence of magnetic field depends on the relation 
between strengths of the fields and the angle between 
them. It is worth to mention the result of acceleration in 
the crossed homogeneous fields. In this case one can find 
the reference frame where the fields are parallel. The 
particle in this reference frame is infinitely accelerated 
along the parallel fields and lose the perpendicular com- 
ponent of the motion. Taking the motion along the field 
and transforming the velocity back to the laboratory ref- 
erence frame, the pitch angle can be expressed as 
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where p = E/B is the ratio of the electric and magnetic 
fields, and 9 is the angle between them. Assuming that 
typically the electric field is smaller than magnetic field, 
we obtain 
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which states that the pitch angle equals the drift velocity 
due to electric drift (in the units of the speed of light). 

In a more general case, a drift due to the electric field 
should appear as well. Then we arrive at a quite general 
conclusion that the radiation of the accelerated particle 
could be considerably different from the curvature radi- 
ation if the electrical drift exceeds the curvature drift 
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The so-called outer gap model of pulsars gives an ex- 
ample of existence of a perpendicular component of the 
electric field. In the acceleration gaps, the electric field 
could not be parallel to the magnetic field everywhere, in 
particular close to the border of the gap the perpendic- 
ular compo nent of the electric field is increased (Cheng 
et al\ |l986| ) . In these regions the radiation can be ditt'er- 
ent from the curvature radiation. 

3. LOCAL TRAJECTORY 

The accurate calculation of the radiation spectrum re- 
quires detailed knowledge of the trajectory of charged 
particle. The latter obtained in the drift approximation 
is not suitable for this purpose since the fast gyrations 
are erased in the course of averaging. Fortunately, under 
the assumption that the gyroradius is small, it is pos- 
sible to obtain a local solution which takes into account 
fast gvrat ions for the inotion in an arbitrary curved mag- 
netic field Trubnikov ( 2000 ) . This approach is equivalent 



to the consideration of particle motion in the magnetic 
field which is constant along binormal and has a con- 
stant curvature and zero torsion. The magnetic lines of 
this field present circles with centres lying on a straight 
line. The magnetic field of this structure is created by the 
current of an infinitely long straig ht wire. It allows us to 



consider the solution obtained in Kelner and Aharonian 
(2012) for such field as a local one. I'his consideration 
is also possible because a small part of the curved mag- 
netic field line can always be approximated by the arc of 
a circle with the radius equal to the curvature radius of 
the line. 

We consider the particle motion in the local coordinate 
system: 

r = sh + xn + yk, (10) 

where h — B/B is the unit vector in the direction of the 
magnetic field, k = h x n is the binormal vector, n — 
rQ(h\7)h is the normal vector wi th rg as the curvature 
r adius . Then, in accordance with Kelner and Aharonian 



( 2012 ) , the local velocity can be written as 
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where 
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/3f "~7' (12) 
Vx = cl3± sinr, Vy = c/Bd — cl3± cost. 



Here is the component of the velocity along magnetic 
field line, /3± is the component of the velocity perpendic- 
ular to the drift trajectory, is the drift velocity defined 
in Eq. (^) (all velocities are in units of c), r = f)t is the 
time in the units of the gyration period 

The expressions in Eq. ([l^ ) describe the motion along 
the magnetic field line with gyration around it and the 
drift in the direction of the binormal to the magnetic 
field line. The solution formally includes also the case of 
strict motion along the drift trajectory without gyration 
{l3± — 0). Such a situation can occur only locally. Since 
the curvature of the magnetic field lines changes, the 
gyration does not disappear completely although it could 
be very small compared to the drift. Eq. is correct if 
the radius of gyration is much smaller than the curvature 
radius. This is equivalent to the conditions f3± <C /3|| and 
/3d ^ /3||- Moreover, it allows us to neglect the drift 
gradient which, otherwise, would lead (in the case of the 
vacuum magnetic field V x B = 0) to the drift velocity 
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The expression for the acceleration in the magnetic 
field 

a=-^(vxB), (14) 
mc7 

and Eq. (|l^) allow us to find the absolute value of the 
acceleration 
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where ag = c^/?|/ro is the acceleration due to curva- 
ture of the magnetic field, ij = (3±/(3d is the ratio be- 
tween the velocity perpendicular to the drift trajectory 
and the drif t velocitv. It is conven i ent to introduce, as 



suggested in Kelner and Aharonian (2012), a parameter 
which shows the difference between the total accelera- 
tion and the acceleration induced by the magnetic field 
curvature 
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Because of the simple relation between the acceleration 
and the curvature radius of the trajectory in the ultra- 
relativistic case a = / Rc, the g-parameter indicates 
also the difference between the trajectory curvature and 
the curvature of the magnetic field lines (compare with 
Eq. (I)) 

Rc = -■ (17) 
9 

As discussed in Sec. 2, we will neglect the difference be- 
tween the curvature of the drift trajectory and the cur- 
vature of the magnetic field line. 

The solution given by Eq. (12) can be applied locally 
if the parameters of the motion such as the energy of the 
particle, the curvature and the strength of the magnetic 
field are changed slowly. In this case, the parameters of 
the solution J7, f3± and /3d change slowly as well. The 
g-parameter varies in the limits |1 — ry| ^ q ^ 1-1-??. 
If 77 ^ 1, then q appears within the band of width 2r/ 
around (7 = 1 . If 77 > 1 , then q varies in the band of width 
2 around q — rj. So over the period of gyration, q could 
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not change greater than 2. Naturally, if the parameters of 
the motion change considerably over the gyration period, 
the behaviour of g-parameter would be different. 

4. RADIATION SPECTRUM 

The spectral power density of the synchrotron radia- 
tion is defined by the i nstantaneous cur vature radius of 
the particle trajectory ( Schwingei 1949). For the curva- 
ture radius given by Eg' ([l7|) the spectral power density 
of radiation of the particle moving in the curved magnetic 
field can be written in the form 



where 
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Here — 3c7'^/(2ro) is the characteristic frequency of 
the curvature radiation and 



(20) 



F{x) K5/3{u)du 



is the emissivity function of synchrotron radiation. The 
g-parameter defined as the ratio between the total accel- 
eration and the acceleration induced by the curvature of 
the magnetic field, can be expressed as 
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where a is the pitch angle, v is the particle velocity, U|| is 
the velocity along the magnetic field. If the energy and 
the magnetic field are changed slowly, the local represen- 
tation given by Eq. (16) can be used. Then the spectrum 
averaged over the period of the gyration is determined by 
the function 



q{ri,T)F 
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In all cases under consideration the form of the spec- 
trum is defined by the same function F{x). The relevant 
parameters change only the position of the maximum 
and the intensity. However, the function F ix) should 
be changed to its quantum analogue (see Eq. |A7| ) if the 
parameter x = -B7 sin a/Bcr is of the order of unity and 
higher, where B cr = 2m'^c^/3eh w 2.94- lO^^ G (jBerestet- 
3kii et al. 1982| ). Note that at such conditions, the energy 



ol the produced photon is close to the energy of the ra- 
diating electron. The electron-positron pair production 
by a gamma-ray photon in the strong magnetic field oc- 
curs when the parameter x 1j where 7 is the photon 
energy in the units of the electron rest mass, and a is 
the angle between the photon and magnetic field. In 
the curved magnetic field the angle between photon and 
magnetic field could become sufficiently large for produc- 
tion of electron-positron pairs. This could lead, provided 
that the optical depth is large, to development of elec- 
tromagnetic cascade and formation of radiation which is 
considerably different from the initial one. 

The substitution of the g-parameter in the general form 
of Eq. ( |2l| ) to Eq. ( p^ results in the standard spectral 
power density of the synchrotron radiation. Thus the 



g-parameter expresses the difference between the con- 
ventional curvature radiation (when the curvature of the 
trajectory is accepted to be equal to the curvature of the 
magnetic field line) and the actual radiation which is the 
small pitch angle synchrotron radiation. One can see that 
there is no well-defined boundary between the curvature 
and the synchrotron radiation. It seems natural to de- 
fine the magnetic bremsstrahlung as curvature radiation 
when the main contribution is introduced by the curva- 
ture of the magnetic field line. This case corresponds 
to r] <^ 1; see Eqs. (^ and ([l^). Then the synchrotron 
radiation occurs when the curvature of the trajectory in- 
duced by the strength of the magnetic field provides the 
main contribution to the radiation. This corresponds to 
(7 « 77 ^ 1. Finally, at ~ 1 both the strength and the 
curvature of the magnetic field play equal role in pro- 
duction of emission which we will call synchro- curvature 
radiation. 

The limits of applicability of Eq. (|8|) defining the en- 
ergy spectra of the synchrotron and the curvature radi- 
ation originate fro m the appr o ach o f its derivation pro- 
posed by Swinger (Schwingei (1949)). The radiation of 
the ultrarelativistic particle is concentrated in a narrow 
cone with the opening angle ^1/7 and is therefore col- 
lected while the angle between the velocity and the di- 
rection of the observation is of the order of the same of 
^ 1/7. The Swinger method is based on expansion of 
the trajectory in the small time interval. During this 
time interval, the entire observable radiation should be 
collected. It means that the approach works if the parti- 
cle velocity changes the direction at an angle larger than 
1/7 while the expansion is valid. The analysis of the lo- 
cal trajectory given by Eq. (|l^ ) in the curved magnetic 
field gives the following limits of applicability: 
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(23) 
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The first condition corresponds to the case of syn- 
chrotron radiation and states that the perpendicular mo- 
tion should be relativistic as it is expected from the 
consideration of radiation in the homogeneous magnetic 
field. The second condition corresponds to the situa- 
tion for the curvature radiation of a particle with a small 
perpendicular momentum. Small gyrations almost do 



not influence on the applicability of Eq. (18) and in the 
limit I3±_ — the condition simply states tr 
should be relativistic. 



aat the motion 



5. NUMERICAL IMPLEMENTATION 

The analytical approach described above allows us to 
study the local properties of the particle trajectory and 
the radiation in the curved magnetic field. To solve the 
problem in the general case, taking into account the en- 
ergy losses of particles, we performed numerical integra- 
tion of the equations of motion. The radiation proper- 
ties have been studied in the dipole magnetic field which 
seems to be a quite good approximation for the strong 
magnetic fields in compact astrophysical objects. The 
dipole magnetic field has two distinct features to be taken 
into account. The first one is the fast decrease of the 
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strength with the distance ~ with a strong impact 
on the radiation intensity. The second one is the signifi- 
cant variation of the curvature with a change of the polar 
angle 9 from the dipole axis cx smO/r, so the radiation 
spectra in the vicinity of the pole and the equator should 
be different. 

For numerical calculations we use the equations of mo- 
tion in the ultrarelativistic limit. In this case the ra- 
diation reaction force is opposite to the velocity which 
changes only its direction. The equation of motion can 
be written in the form 



mc-(7/3) = e(/3 x B) - |/|/3, 



(24) 



where (3 is the velocity in units of c with \(3\ w 1, J5 is 
t he magnetic field, and / i s the radiation reaction force 
( Landau and Lifshitz 1975). Taking the scalar product 
of Eg. (|2^) with velocity (3 we obtain the differential 
equationTor the energy losses. This equation allows can- 
cellation of the radiation reaction force with the Lorenz 
factor time derivative in Eq. (p^). Finally, the equa- 
tion of motion has the same form as for the considera- 
tion without energy losses, where the energy enters as a 
parameter. For the sake of convenience of the numeri- 
cal treatment and comprehension of the structure of the 
system of equations, these equations are written in the 
dimensionless form: 



dr 



= wi/3, 

= e(/3 X b) , 

= M2(/3 X bf. 



(25) 
(26) 
(27) 



The system of equations depends on two dimensionless 
parameters 
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where 70 is the initial Lorenz factor, i?o is the charac- 
teristic distance to the radiating region from the dipole, 
Bq = B^, [R^/ Rof' is the characteristic magnetic field 
with i?* and being the source radius and the mag- 
netic field at the pole of the dipole, to is a particle mass, 
c is the speed of light. Here we have introduced the fol- 
lowing dimensionless variables: r' is the coordinate in the 
character units of length Rq, (3 in the velocity of the par- 
ticle in the units of c, ^ = 70 /7 is the ratio of the initial 
and current value of the Lorenz factor, r = teBo/mcyo 
is the characteristic time in the units of the initial gyra- 
tion period, b is the dimensionless dipole magnetic field 
which is expressed as 
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(29) 



where fi is the unit vector in the direction of the dipole 
axes, n = r'/r' is the unit vector to the particle position. 

It should be noted that depending on the specific condi- 
tions characterizing an astrophysical source, the param- 
eters ui and U2 may differ by many orders of magnitude. 
For instance, for typical parameters of the so-called po- 
lar cap model of pulsar magnetosphere Bq = 10^^ G, 



i?o — 10^ cm, and 70 = 10^, we have ui ~ 1.7 • 10"'' and 
U2 ~ 1.1 • 10^2. Thus the problem is non-stiff, therefore 
for integration of this system of differential equations the 
implicit Rosenbrock method has been used. 

The calculations of the trajectory have been carried out 
for different initial conditions. The initial position deter- 
mined by the radius Rq and the polar angle 6*0 relative to 
the magnetic dipole axis defines the typical environment 
parameters for the models under consideration. The ra- 
diation spectrum has been studied for different initial 
pitch angles ao and initial Lorenz factors 70. The de- 
tailed knowledge of the trajectory and the energy allows 
us to find g-parameter from Eq. (pT|) and the radiation 
spectrum given by Eq. ( |l8|) at any moment of time. The 
spectra integrated over the time, the so called cumulative 
spectrum, have been obtained for any initial condition 
under consideration. 

Finally, in the case of a very strong magnetic field 
and/or very large Lorenz factor, the particle may radiate 
in the quantum regime. More specifically, when the pa- 
rameter X = B"f sina/ Bcr becomes of the order of unity 
or larger, Eq. (^ should be replaced by its quantum 
analogue (see Appendix) 



dr 



M2(/3 X 6)2 



(1+U3((/3X b)/e)2/3 



(30) 



where U3 = 1.07- 10-9(So7o)^/^ and B^ 
2.94-10" G. 



6. ASTROPHYSICAL IMPLICATIONS 

In this section we explore possible realizations of the 
synchrotron and curvature regimes, as well as transi- 
tions between these two modes of radiation (the synchro- 
curvature regime) in the context of three specific astro- 
physical scenarios. Namely below we discuss the radi- 
ation of electrons and positrons in the pulsar magneto- 
sphere for the outer gap and polar cap models, and the 
radiation of protons at their acceleration in the vicinity 
of a supermassive black hole. 

6.1. Outer Gap 

The high energy gamma radiation from pulsars is 
widely believed to be produced in the ou t er gap of the 



pulsar magnetosphe re ( |Cheng et al. 1986; Takata et al. 



2004; Hirotani) 2005). Here we present the results of our 



calculations of radiation of electrons (positrons) in the 
dipole magnetic field at the location of the outer gap. 
The position of the outer gap was assumed along the last 
open field line which is is determined by the inner bound- 
ary located at the null surface where the Goldreich- Julian 
charge density is zero, and by the outer boundary taken 
on the surface of the light cylinder. We adopt the param- 
eters of the Crab pulsar: the radius of the star R^, = 10^ 
cm, the rotation period P = 33.5 ms, and the magnetic 
field at the pole B^ = 10^^ G. We consider non-aligned 
pulsar with an angle of 7r/4 between the rotational axis 
and the magnetic dipole axis. The particle initial posi- 
tion is set at the distance rinu = 0.5Ric from the null 
surface along the last open field line, where Ric is the 
radius of the light cylinder. During the numerical cal- 
culations, the particle is followed up to the intersection 
with the light cylinder. For the Crab parameters these 
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conditions correspond to 9q — 51.5°, Rq = 1.2 • 10* cm, 
Bo = 5.6 • 10^ G, and Ric = 1.6 • 10* cm, where 6*0 is the 
polar angle relative to the magnetic dipole axis. 

The spectra and the radiation regimes have been stud- 
ied for different initial directions and Lorenz factors of 
electrons. The resulting spectra are shown in Fig. |l|. The 
complementary plots in Figs. |, ||, and ^ demonstrate the 
time-evolution of the value of g-parameter (blue lines, 
right scale) and the Lorenz factor normalized the initial 
Lorenz factor. The local behaviour of the q-parameter 
agrees well with Eq. (|6|) obtained from the local solu- 
tion. As discussed above, the oscillations caused by the 
particle gyration are within the bandwidth of < 2. The 
complementary plots allow us to observe simultaneously 
the energy lose rate and the regime of the radiation pre- 
sented by the value of the q-parameter. If the energy loss 
rate is low, the rarefied output of the points produces in- 
termittent curves for the q-parameter. 

It is interesting to examine the statement that the par- 
ticle moving along drift trajectory yields minimum inten- 
sity of radiation and produces less energetic photons. To 
do this, the initial velocity is deflected by the angle 
from the magnetic field line towards the binormal vector. 
The corresponding cumulative (integrated along trajec- 
tory) energy spectra of radiation are shown in Fig. |^. 
The complementary plots are shown in Fig. |^. We can 
see that indeed for the initial Lorenz factors 70 = 10^ 
and 7o — 10*" the radiation is less energetic compared to 
other cases and the rates of energy losses are minimal as 
well (compare with corresponding plots in Figs. ^ and|^). 

The situation however is different for larger initial 
Lorentz factors; see the curves corresponding to 70 = 
5 • lO"^ and 70 — 10*. It is seen that for the initial di- 
rection along the drift trajectory the radiation spectra 
extend to higher energies than in the case of the initial 
direction along the magnetic field. This can be explained 
by the very intensive energy losses occurring even before 
the particle has made the first gyration (the first oscil- 
lation of g-parameter). During this time q ~ 1 for the 
initial direction along the drift trajectory and q ss for 
the initial direction along the magnetic field. Eqs. (|l8|)- 
(19r) show that larger values of q give a more energetic 
radiation. Correspondingly, more energetic radiation is 
produced in the case of the initial direction along the 
drift trajectory. However, after many gyrations the en- 
ergy losses in the case of the initial direction along drift 
trajectory become, as expected, less intensive compared 
to the case of the initial direction along magnetic field 
line. 

We call the attention of the reader to the regimes of 
the radiation demonstrated by curves in Fig. 0. For 
70 = 10^ and 70 = 10^, the q-parameter equals unity 
indicating that the radiation proceeds in the curvature 
regime. However g = 1 is not exact as it is demonstrated 
for 7o — 10^ where the scale for g-parameter slightly 
oscillates around unity. These small oscillations corre- 
spond to the fine gyration around the drift trajectory. 
As discussed above (see Eq. |2^) the presence of such fine 
perpendicular motions does not infiuence on the applica- 
bility limits. 

In the cases 70 = 5 • 10^ and 70 = 10* the q-parameter 
has more complex behaviour. The increase at the be- 
ginning is defined mostly by the fast energy losses. The 
decrease is defined by the combination of several fac- 



tors such as the reduction of the magnetic field strength 
and the change of its curvature. The increase of q- 
parameter indicates that the radiation occurs in the 
synchro-curvature or the synchrotron regimes when the 
radiation due to curvature of the magnetic field line is 
less important (for 70 = 5 • 10^ or simply negligible (for 
7o = 10*). According to Eq. (p^, the energy of the radia- 
tion maximum scales as ~ J^q- The g-parameter reaches 
the maximum when the considerable amount of energy 
has been lost. Therefore, in spite of large q, the peak of 
radiation shifts towards low energies and does not affect 
the cumulative spectrum. The interesting feature can be 
seen at first moments when the energy oscillates with q- 
parameter and the minimums of g-parameter correspond 
to fiatter parts of the Lorenz factor evolution curve. 

The radiation spectra of electrons launched along the 
magnetic field line are slightly more energetic except for 
7o = 5 • 10^ and 70 — 10*. The reason is the same as 
discussed above. Initially, the radiation for 70 = 10^ and 
7o — 10^ is in transition regime, although for 70 — 10^ 
the most of the energy is lost in the curvature regime 
which occurs fast. Therefore the spectrum in this case 
almost coincides with the spectrum for the initial direc- 
tion along the drift trajectory. For 70 = 10^ the most 
of the energy is lost in the transition regime, thus the 
spectrum is shifted to higher energies by g « 2. 

For illustration of the effect related to initial pitch an- 
gles larger than , we show the case of initial direction 
deflected at the angle lO/?^) from the magnetic field line 
towards the direction opposite to the normal vector. The 
corresponding spectra indicated in Fig. |l| by dash-dotted 
lines are shifted towards higher energies. Although the 
g-parameter (Fig. ^) reaches large values, the most of 
the energy is lost at initial stages at g « 10. Thus the 
spectra are shifted by g « 10 compared to the curvature 
radiation spectra. 

We should note that the energy spectra of radiation 
produced in all regimes contain an exponential cut-off 
at highest energies similar to the spectrum of the small- 
angle synchrotron radiation. However, since the radia- 
tion spectrum is very sensitive to the pitch-angle, a pop- 
ulation of electrons with similar energies but different 
angles can result in a superposition spectrum with a less 
abrupt cut-off. The condition for realization of such a 
spectrum is that the distribution of electrons over pitch- 
angles around zero angle should be wider than f3o- 

6.2. Polar Cap 

In the polar cap model electrons radiate in the region 
located close to the surface of the neutron star where 
the magnetic held is much stronger than in the outer 
gap model, approaching to B Ri 10^^ G. This results in 
much faster damping of the perpendicular component of 
motion. The very small drift velocity /3d implies that 
the drift trajectory and the magnetic field line almost 
coincide, thus even a small deflection from the magnetic 
field line produces radiation quite different from the cur- 
vature radiation. However, the transition to the curva- 
ture radiation regime occurs very fast. In the curvature 
regime, electrons radiate more energetic photons than in 
the outer gap model. But the curvature of the magnetic 
field lines in the polar cap model with ~ 1° is only an 
order of magnitude larger than in the outer gap since the 
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Figure 1. The cumulative (integrated along trajectory) radiation 
spectra of electrons calculated for the outer gap model of the pulsar 
magnetosphere. The curves are obtained for different initial Lorenz 
factors of electrons 7 = 10^,10^,5 X 10^, 10®, and for different 
initial directions relative to the magnetic field lines: along the drift 
trajectory (solid lines), along the magnetic field line (dashed lines), 
and at pitch angle a = lO/S^i in the meridional plane opposite to 
the normal vector of magnetic field lines (dashed-dotted lines). 
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Figure 2. Time evolution of the g-parameter and the electijon 
Lorenz factor in the outer gap model (complementary to Fig. 111). 
Four panels correspond to the initial Lorenz factor of electrons 
7 = lO'^, 10^, 5 X 10 , 10^ and their initial direction along the drift 
trajectory. 

curvature of the dipole magnetic field scales as ~ sinO/r. 
Correspondingly, the maximum energy of curvature ra- 
diation at the polar cap is only an order of magnitude 
higher than in the outer gap. 

The energy spectra of radiation calculated for the po- 
lar cap model are shown in Fig. ^ The complementary 
plots for the evolution of the g-parameter and the elec- 
tron Lorenz factor are presented in Figs. ^, ^ and |^. 
The initial position of the particle is i?o = 10 cm, and 
00 = 3°. The spectra indicated by solid lines correspond 
to the case when the initial direction of particle is along 
magnetic field line. In this case the radiation is in the 
curvature regime. But at the initial stage, the radia- 
tion proceeds in a very fast synchro-curvature regime. 
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Figure 3. The same as in Fig.H, but for the initial direction of 
electrons along the magnetic fielctline. 





Figure 4. The same as in Fig.^, but for the initial direction of 
electrons at an angle 10/3 jj. 

but abruptly turns to the regime with 5 « 1 and fast 
oscillations caused by fine gyrations. 

The abrupt change of the regimes leads to an interest- 
ing feature in the cumulative spectra for the case of an 
initial pitch angle uq — I/70 (dashed lines). Because 
of small changes in energy and fast change of the q- 
parameter, the energy spectra consist of two peaks. The 
peak at higher energies is produced by synchrotron ra- 
diation (g ^ 1, see Fig. |^), while the lower energy peak 
corresponds to curvature radiation regime. The double- 
peak structure disappears for large initial pitch angles. 
For example, for the pitch angle IOO/70 the transition to 
the curvature regime is very fast, and the electron enters 
into this regime with dramatically reduced Lorentz fac- 
tor. Thus the peak of the curvature radiation not only 
is shifted to smaller energies, but also is too weak to be 
seen in the cumulative spectrum^. 

In very strong magnetic fields, namely when the pa- 

^ Note that the pitch angle IOO/70 which we treat as 'large', still 
is extremely small, ~ 2(7o/10^)~^arcsec. 
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Figure 5. The cumulative (integrated along trajectory) radiation 
spectra of electrons calculated for the outer gap model of the pulsar 
magnetosphere. The curves are obtained for different initial Lorenz 
factors of electrons 7 = lO'', 10^, 10®, and for different initial direc- 
tions relative to the magnetic field lines: along the magnetic field 
line (dashed lines), and for two pitch angles I/70 (dashed lines) 
and 100/70 (dashed-dotted lines). 
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Figure 6. Time evolution of the g-parameter and the electmn 
Lorenz factor in the outer gap model (complementary to Fig. H). 
Three panels correspond to the initial Lorenz factor of electrons 
7 = 10^,10^^,10* and their initial direction along the magnetic 
field line. 



rameter x — B'f sin a /Bcr > 1, the radiation is produced 
in the quantum regime. Let's assume that the initial 
pitch angle is inverse proportional to the initial Lorentz 
factor, ao — a/70. This makes the condition of radiation 
in the quantum regime independent on 79: 



X 



B 



-7sm ao 



1 • W^G 
2.94- IQI^g' 



3.4- 10" 



(31) 



Thus, at the initial pitch angle with a > 30, the elec- 
trons radiate in quantum regime. Dash-dotted lines in 
Fig. 1^ present radiation spectra for the initial pitch- 
angles Qfo = 100/70. Note that in the quantum regime 
almost the entire energy of the parent electron is trans- 
ferred to the radiated photon. Thus we should expect 
abrupt cutoff in the radiation spectra. This effect is 
clearly seen in Fig. || (dot-dashed curves corresponding 
to the initial pitch angle ao = IOO/70). The gamma rays 
produced in the quantum regime are sufficiently ener- 
getic to be absorbed in the magnetic field through the 
e"'"e~ pair production. This will lead to the development 
of an electromagnetic cascade in the magnetic field. The 
spectrum of cascade gamma-rays that escape the pulsar 
magnetosphere will be quite different from the spectra 
shown in Fig. ^. 

6.3. Protons in the Black Hole magnetosphere 

The acceleration of the protons to ultrahigh energies 
in the potential gap of the spinning supermassive black 
hole should be accompanied by curvature radia tion in 
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Figure 7. The same as in Fig.g, but for the initial direction of 
electrons at pitch angle 1/70- 

The larger (compared to electron) mass of proton leads 
to the larger drift velocity Therefore the spectra 

and the radiation regimes in the case of protons are less 
sensitive to the initial direction of motion. The radiation 
of protons deviates from the curvature radiation when 
the perpendicular component of proton's Lorentz factor 
exceed 



7-L = 7/3d 



9 '? 
rripC J 

eBro 



= 2.7-10 



Vioio/ 



(32) 



( p002[) ). Here we briefly 



the r nagne ti c field which thre a ds the black hole ( |Levin-| This is larger by the factor of mp/rrie « 2 • 10^ com- 
pared to the same condition for electrons. It is inter- 
esting to note that for the same Lorenz factor, in the 
synchrotron regime protons radiates much weaker than 
electrons, whereas in the curvature regime they radiate 
equally. 

The radiation spectra of protons with different initial 
pitch angles are presented in Fig. ^. The smallest angle 
ao is close to /3d, and the proton radiates in the synchro- 



Aharonian et al 
examine the radiation regimes and the gamma ray spec 
tra of accelerated protons produced in this model. As 
before, we consider the radiation in the dipole magnetic 
field adopting the following (typical for a SMBH) param- 
eters -B* = 10^ G and i?* = 10^^ cm. The initial position 
of the particle is set at the polar angle 6*0 = 5° relative 
to the magnetic dipole axis. 



10"^° 



Figure 8. The same as in Fig.|^, but for the initial direction of 
electrons at pitch angle 100/70- 

curvature regime (solid lines). Because of high energy 
the gyration period is large, and in the case of 70 = 10^'^ 
the particle makes only several gyrations before escaping 
the region with high magnetic field (see Fig. |l^). After 
that, the proton radiates with a very low rate, and the q- 
parameter approaches zero. In the case of 70 = 10®, the 
proton gyrates more frequently (see Fig. ^ and there is 
seen more pronounced tendency of q approaching to zero. 
The other curves correspond to the initial pitch angles 
lOao and lOOao. The spectra are shifted towards higher 
energies by the same factors of q ~ 10 and q w 100. Ac- 
cordingly, for the Lorenz factor 70 = 10^^ which is not 
shown in figures, the proton will radiate predominantly 
in the curvature regime and the spectrum will be shifted 
by a factor of (lO^/lOi'')^ = 1000 relative to the cur- 
vature spectrum of the proton with initial Lorenz factor 
7o = 10^". Despite the small energy losses relative to 
the initial energy, the amount of the radiated energy is 
quite large. As in the case of the pulsar polar cap model, 
for the chosen parameters the curvature radius of the 
magnetic field lines is larger by an order of magnitude 
compared to the gravitational radius of the black hole 
which usually is taken for evaluation of curvature radius. 
It yields smaller energy losses and increase the maximum 
Lorenz factor of acceleration compared to the case when 
the gravitational radius of the black hole is used as a cur- 
vature radius. To be more specific, the maximum Lorenz 
factor of a particle in radiative-loss limited regime scales 



with c u rvature radius R r as Rd R* (Aharonian et al 
( p002| ), |Levinson| ( |2000| )). 



7. SUMMARY 

The radiation of relativistic particles in a strong curved 
magnetic field is of a great astrophysical interest, in 
particular in the context of magnetospheric gamma-ray 
emission of rotation-powered pulsars. In these objects 
with very strong magnetic fields, the radiation proceeds 
in quite complex regimes, so it cannot be reduced merely 
to the consideration of the nominal synchrotron and cur- 
vature channels. For proper understanding of these ra- 
diation regimes and the transitions between them, the 
accurate treatment of the particle trajectory is a key is- 
sue. It can be done by solving the equations of motion 
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Figure 9. The cumulative (integrated along trajectory) radiation 
spectra of protons in the magnetic field of rotating supermassive 
black hole. The chosen parameters are described in the text. The 
spectra correspond to two different Lorentz factors and three dif- 
ferent initial pitch angles relative to the magnetic field lines: 10~* 
(solid), 10"^ (dashed), 10"^ (dash-dotted) radian for 70 = lO'' 
(red) and 10"^ (solid), 10"^ (dashed), 10"^ (dash-dotted) radian 
for 70 = 10^'' (blue). 
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Figure 10. Time-evolution of the Lorenz factor and the q- 
parameter in the magnetic field of a supermassive black hole with 
the initial Lorenz factor of protons 70 = 10^ and their initial pitch 
angles a = 10""', 10"^, and 10'^. 

taking into account energy losses of charged particles. In 
this work we followed this approach to explore the ra- 
diation features of ultrarelativistic particles in a dipole 
field which is a good approximation for the magnetic field 
structure in compact astrophysical objects. The accu- 
rate numerical solutions for the test particle trajectory 
allowed us to trace the radiation regimes and calculate 
self-consistently the radiation spectrum without any a 
priori assumption concerning the radiation regime. 

In this paper we demonstrate that even small deflec- 
tions of the initial particle motion from the magnetic 
field lines may result in a radiation spectrum quite dif- 
ferent from the spectrum the curvature radiation. For 
any initial pitch angle, the particle tends, while losing its 
energy, to turn to the curvature radiation regime and to 
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Figure 11. Time-evolution of the Lorenz factor and the q- 
parameter in the magnetic field of a supermassive black hole with 
the initial Lorenz factor of protons 70 = 10"'^'' and their initial pitch 
angles a = 10~^, lO""^, and 10~-'^. 



move strictly along drift trajectory (although never can 
achieve the latter). In principle, the fast increase of cur- 
vature of the magnetic field lines can turn the regimes in 
inverse order. 

In different environments (or different positions rela- 
tive to magnetic dipole) the transition to the curvature 
radiation regime proceeds with different rate. While in 



the polar cap model the transition occurs almost in- 
stantly, in the outer gap the transition could last so long 
that the particle would not turn to the curvature regime 
while passing the gap. The spectrum of the radiation be- 
comes very different from the spectrum of the curvature 
radiation if the pitch angle of the particle is greater than 
drift speed /3d — vd/c^ and quite similar if it is smaller 
than Pd. For the typical parameters of the polar cap 
model, = 2.2-10-9 (7/1 0^). This implies that the 
even tiny deflection from the magnetic field line leads 
to the spectrum different from the curvature radiation 
spectrum. 

Significant deviations of the radiation spectra from the 
nominal curvature radiation spectrum is expected also 
in the outer gap model. The particles do not move 
along magnetic field lines but gyrate around drift tra- 
jectory. This results in a different (more "energetic") 
spectrum from the conventional curvature radiation spec- 
trum. In the outer gap model, electrons with initial di- 
rection along magnetic field line start to radiate in the 
synchro-curvature regime, when the both the curvature 
of the magnetic field and its strength play equally impor- 
tant role. The effect is quite significant at large Lorentz 
factors of electrons, 7 > 10^. 

Finally, we demonstrate the strong impact of the initial 
pitch angle on the radiation spectrum in the scenario of 
acceleration and motion of ultrahigh energy protons in 
the magnetosphere of a supermassive black hole. 

A. P. is a fellow of the International Max Planck Re- 
search School for Astronomy and Cosmic Physics at the 
University of Heidelberg (IMPRS-HD). 



APPENDIX 

ENERGY LOSSES IN THE QUANTUM REGIME 
In a very strong magnetic field, the ultrarelativistic electrons can radiate in the quantum regime, provided that 



B 

X = — 7sma > 1, 

JDc.r 



(Al) 



where B^r = 
et a/.||l982| ) 



3e?i 



2.94 • 10 G. The energy lose rate can be written in the form (Bayer et al. 1973; Berestetskii 



dE 



where 



and 



dt 



H{x) 



H{T,x)dT , 



H{t, x) = X [(1 - r)Fix) + XT^K2/3ix)] , x 



X(l-r)' 



(A2) 
(A3) 
(A4) 



where -^^2/3 (2;) is the modified Bessel function of the order 2/3, F{x) is the emissivity function of the synchrotron 
radiation (see Eq. |2^), t = e/E, where e is the energy of the radiated photon, E is the energy of the radiating particle. 

For calculations it is convenient to express H{x) in Eq. ( |A2| ) in a simple approximate analytical form. Using asymp- 
totics of this function 



i7(x)«^2^/3rmx 

243 1 3'-^ 



(A5) 
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we have found the foUowing approximation 

87r\/3 



H{X) 



27 



^ , 3 (2x)^/3 

r(|), 



X 1 + 



0.52Vx(l + 3Vx-3.2x) 
l + 0.3Vx+17x+llx2 



(A6) 



The first two terms (before the sign x ) of Eq. ( A6) give right asymptotics at x ^ 1 and x ^ 1 and provide an accuracy 
better than 10% for other values of x, whereas the inclusion of the last term in the brackets makes the accuracy better 
than 0.1% for any x- 
The spectrum of the radiation in the quantum regime is expressed as 



Fq{x, r) = (1 - t)F{x) + t'xK2/3{x), x= — 



E 



(AT) 



where = '^^^2mc ° '^^ ^^'^ characteristic energy of the emitted photon. To use this function in Eq. ([l^), the ec 
should be changed to /la;,. An analytical approxi mation of this function can be obtained using the approximation for 
emissivity function of the synchrotron radiation ([Aharonian et aJ. 2010) 



F{x) w 2.15a;i/3(l + SMx)^/ 



i/gl + 0.884a;2/3 + 0.4715 



.4/3 



1 + 1.64x2/3 + 0.974x4/3 



and 



xK2/3{x) w 1.075x^/^(1 + 3.72x)^/ 



i/gl + 1.58x^/3 + 3.97x4/3 _^ 



1 + 1.53x2/3 + 4.25x4/3 

Both approximations provide an accuracy better than 0.2% for any value of the argument x. 
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